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Abstract 

We present new groundbased 3 fim spectra of 14 young stellar objects with H 2 ice 
absorption bands. The broad absorption feature at 3.47 fim was detected toward all objects 
and its optical depth is correlated with the optical depth of H 2 ice, strengthening an 
earlier finding. The broad absorption feature at 3.25 fim was detected toward two more 
sources and an upper limit is given for a third source. The optical depths of the 3.25 
^m feature obtained to date are better correlated with the optical depth of the refractory 
silicate dust than with that of H2O ice. If this trend is confirmed, this would support 
our proposed identification of the feature as the C-H stretch of aromatic hydrocarbons at 
low temperature. An absorption feature at 3.53 fim due to solid methanol was detected 
for the first time toward MonR2/IRS2, as well as toward W33A and GL 2136. The 
wavelengths of the CH 3 OH features toward W33A, GL 2136, and NGC7538/IRS9 can be fit 
by CH 3 OH-rich ices, while the wavelength of the feature toward MonR2/IRS2 suggests an 
H 2 0-rich ice environment. Solid methanol abundances toward GL 2136, NGC7538/IRS9, 
and MonR2/IRS2 are 3-5 % relative to H 2 ice. There is an additional narrow absorption 
feature near 3.47 /im toward W33A. For the object W51/IRS2, spatially resolved spectra 
from 2 to 4 /im indicate that the H 2 ice is located predominantly in front of the eastern 
component and that the H 2 ice extinction is much deeper than previously estimated. 
For the object RNO 91, spectra from 2 to 4 /im reveal stellar (or circumstellar) CO gas 
absorption and deeper H 2 ice extinction than previously estimated. 
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1. Introduction 



Several distinct spectral absorption features have been detected in the long-wavelength 
wing of the interstellar 3.1 fim H 2 ice absorption band seen toward molecular cloud 
sources. As the wing extends through the 3.2-3.6 //m region characteristic of C-H stretch 
vibrations, some or all of these features could be due to organic species in the solid state. 
Spectra of embedded young stellar objects in this region have so far revealed a broad 3.47 
Hm feature and a broad 3.25 /xm feature, neither of which is securely identified yet, and 
a feature at 3.53 due to solid CH3OH. There is also a reported detection of a narrow 
feature near 3.47 fj,m, attributed to solid H 2 CO, toward GL 2136 (Schutte et al. 1996a). 
The species responsible for the extended long-wavelength wing itself is still not yet certain. 

The broad 3.47 ^m feature (FWHM^ 0.1 ^m) was first noticed in four objects by 
Allamandola et al. (1992). They suggested that the feature might be due to the C-H stretch 
absorption of solo hydrogens attached to sp 3 bonded carbon clusters, the "diamond" -like 
form of carbon. The feature was present in every molecular cloud source looked at by 
Brooke, Sellgren, & Smith (1996), with an optical depth that appeared to be correlated 
with the depth of the H 2 ice band, as opposed to that of the silicates. The feature has also 
been seen toward Elias 16 (Chiar, Adamson, & Whittet 1996), a star behind the Taurus 
cloud, so the presence of the feature does not appear to require any heating or radiation 
processing by an embedded protostar. Graham (1998) detected the feature in two sources 
in the RCrA association. Under the assumption that the feature is due to C-H bonds, 
Brooke, Sellgren, & Smith (1996) interpreted the correlation with H 2 ice as indicating that 
both C-H bonds and H 2 ice form in step on molecular cloud dust by hydrogen addition 
reactions. However, they noted that other identifications of the feature are also possible. 

The broad 3.25 jum feature was first noted by Sellgren, Smith, & Brooke (1994) 
in a spectrum of MonR2/IRS3, and later confirmed with better spectra (Sellgren et al. 
1995). They speculated that the feature might be due to aromatic hydrocarbons at low 
temperature. Two more detections were reported (Brooke, Sellgren, & Smith 1996), but 
there were insufficient data to determine any correlations with H 2 ice or silicates. The 
feature falls in a wavelength region where there are strong H 2 and CH 4 absorption lines 
in the terrestrial atmosphere. It also falls on the steeply sloping edge of the H 2 ice band. 
Thus reliable detections of the feature are extremely difficult. 

Solid methanol has now been seen toward several embedded sources in molecular 
clouds. Solid CH3OH was first identified in the interstellar medium by Grim et al. (1991). 
They observed an absorption feature at 3.53 ^m toward W33A, in good agreement with the 



-4- 



u 3 symmetric CH 3 stretch mode in solid CH 3 OH. This feature was subsequently detected 
toward two other sources: NGC7538/IRS9 (Allamandola et al. 1992) and GL 2136 (Schutte 
et al. 1996a), and has recently been detected toward another, RAFGL7009S (Dartois 
et al. 1998). Additional features which can be attributed to solid CH3OH have been 
detected toward GL 2136 at 8.9 and 9.7 /im (Skinner et al. 1992) and toward W33A and 
RAFGL7009S at 3.85 and 3.94 ^m (Geballe 1991; Allamandola et al 1992; Dartois et al 
1998). Also, part of the broad 6.8 /im feature seen toward many objects is most likely due 
to solid CH3OH (Schutte et al. 1996b and references therein). The wavelengths of peak 
absorption have been used to infer that solid CH 3 OH is found in ice mixtures which are not 
dominated by H 2 ice (Skinner et al. 1992; Allamandola et al. 1992; Schutte et al. 1996a; 
Dartois et al 1998). 

As infrared instrumentation continues to improve, it has become possible to detect 
weak spectral features in an increasing number of young stellar objects. New 3 fim spectra 
of 14 objects with H 2 ice bands were obtained in 1995 and 1996 with two 256 x 256 InSb 
array spectrometers to determine the prevalence and abundance of organics in molecular 
cloud dust. These were objects with right ascensions 16 h < a < 21 h , with the exception of 
MonR2/IRS2. MonR2/IRS2 was re-observed to obtain a spectrum of better quality than 
that shown in Brooke, Sellgren, & Smith (1996), which suffered from poor cancellation of 
telluric features. 



2. Observations 



Spectra were obtained at the Kitt Peak National Observatory 2.1-m telescope on 12 
January 1996 and 4,6 June 1996 and at the United Kingdom Infrared Telescope (UKIRT) 
on 6 September 1995, 12-13 June and 16 July 1996. All dates are UT. Table 1 is a log of 
the observations. The spectra are shown in Figs. 1-4. 

At Kitt Peak, the infrared cryogenic spectrometer CRSP was used. The pixels were 
0.61" across. The slit was 1.6" wide and oriented N-S for MonR2/IRS2 and E-W for the 
others. Background subtraction for MonR2/IRS2 was done by nodding the telescope 10" 
S along the slit. For the other objects, the telescope was nodded 160" E (W for Serpens 
SVS 20). The spectra were sampled every 1/2 of a resolution element. Observations of 
MonR2/IRS2 and Serpens SVS 20 in the 3.5 /im region were taken with the 300 lines/mm 
grating (Grating 1), which provided a spectral resolution A/AA « 1300 at 3.5 fim. The 
spectrum of SVS20 was smoothed to an effective resolution of 800 to improve the 



- 5 - 



signal-to-noise. Observations of W51/IRS2 and RNO 91 were taken with the 150 1/mm 
grating (Grating 3), which provided A/AA ~ 530 in the 3.4-3.9 jum region and ~ 700 in the 
2.2 /im region. The spectrum of RNO 91 was smoothed to an effective resolution of ~ 200 
in the 3.6-3.9 /xm region and R3 400 in the 2.2 /im region. 

At UKIRT, the facility spectrometer CGS4 was used. The pixels were 1.2" across. 
A 1.2"-wide slit was oriented E-W. The sources were nodded ~ 12" along the slit for 
background subtraction. The spectra were sampled every 1/3 of a resolution element. 
Spectra in September 1995 were obtained with the 75 1/mm grating which provided a 
spectral resolution of ~ 1200. The rest of the spectra were obtained with the 150 1/mm 
grating which provided a spectral resolution of ~ 2400. They were smoothed to an effective 
resolution of R3 1300. 

For those sources observed at 2.2 /im, the telescope was first positioned to the 3.5 
jum peak, since surrounding nebulosity could cause the peak position to be a function of 
wavelength. 

The stars used for atmospheric correction were also used for flux calibration. The stars 
were assumed to be blackbodies. The assumed blackbody temperatures based on spectral 
type are given in Table 1. The absolute flux calibration may be uncertain by 20 %. Where 
possible, Kitt Peak spectra were scaled to the flux level of the UKIRT spectra, which we 
consider more reliable. The scaling factors were less than 20 %. 

Hydrogen absorption lines in the standard stars were not removed. This means that 
apparent emission features at 3.04 /im (HI Pfe) and 3.30 /im (HI Pf<5) in the data for 
Elias 29, GL 2136, SVS 20, RCrA/IRSl, and RCrA/IRS2 may be due all or in part to 
corresponding absorption features in the standards. In the case of RNO 91, the standard 
star used to correct the K(2.2 /im) and 3.7 /im regions observed at Kitt Peak was an A2V 
star with strong hydrogen absorptions at Br7 (2.17 /im) and Pf7 (3.74 /im), so these regions 
were deleted. The emission lines in the data of W51/IRS2(E) are predominantly intrinsic 
to the source. 

Some data points in the strongest telluric absorption features were not well-corrected 
and were dropped. Other gaps in the spectra are due to occasional bad pixels. Error bars 
are omitted in the flux plots for clarity, but representative error bars are shown in the plots 
of optical depth that follow. 



3. 2-4 /im Spectra 
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Two young stellar objects were observed from 2 to 4 y«m in order to obtain improved 
estimates of the H 2 ice band optical depth. These spectra are shown in Fig. 4 and 
discussed in this section. Our goal was not an in-depth study of the 2 //m spectra of these 
objects, but some interesting aspects of the spectra are noted. 



3.1. W51/IRS2 



W51/IRS2 is an infrared source that contains an embedded HII region. It has a 3 fim 
H 2 ice absorption band previously estimated to have an optical depth r « 0.8 (Soifer, 
Russell, & Merrill 1976; Joyce & Simon 1982). These ice band observations were taken 
with aperture diameters of 17" and 11", respectively. The source was shown to be double 
in the mid-infrared by Genzel et al. (1982). The eastern source, which is brighter over the 
2.2 to 20 yum range, has a 10 y«m silicate absorption feature, but the western source ~ 5" 
arcsec away has no detectable silicate feature. Goldader and Wynn- Williams (1994) showed 
that W51/IRS2(E) has a 3 y«m ice absorption band but that W51/IRS2(W) does not. Our 
spectra show that the ice band toward W51/IRS2(E) is deeper than the previous estimate. 

The 2-4 ^m spectrum of W51/IRS2(E) is shown in Fig. 4a. The measured fluxes 
for W51/IRS2(E) are within 20 % of those of Goldader & Wynn- Williams (1994). The 
slits were oriented E-W and centered on W51/IRS2(E). With this orientation, we also 
obtained a spectrum of a bright region about 2" S of the position given for W51/IRS2(W) 
by Goldader & Wynn- Williams. Emission lines of HI, Hel and H 2 were detected toward 
both regions. The detected lines are indicated in the figure for W51/IRS2(E) (Fig. 4a). 
The equivalent width of the HI Br 7 line (corrected for absorption in the standard star) 
in the region 2" S of IRS2(W), W A (Br 7 )^ 540 A, is very close to the theoretical value for 
emission by thermal electrons at T=10,000 K in an HII region 590 A, Wynn- Williams 
et al. 1978), but in W51/IRS2(E) it is lower (W A (Br7)^ 140 A), indicating that the 2 
^m continuum of W51/IRS2(E) contains starlight, scattered light and/or dust thermal 
emission. This is in agreement with the conclusion of Goldader & Wynn- Williams (1994) 
from spectrophotometric imaging. Line ratios are given in the caption to Fig. 4a. 

W51/IRS2(E) exhibits a strong H 2 ice band. There was no evidence for any H 2 ice 
absorption in the western region. The continuum chosen for the W51/IRS2(E) ice band is 
shown in Figure 4a. The other curve is a scaled blackbody to indicate the approximate color 
temperature of 450 K. The local continuum adopted for the 3.47 yum feature is also shown. 
The H 2 ice is partially crystalline as indicated by the sharp absorption peak at 3.1 /im and 
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the inflections at 2.97 and 3.20 pm. The latter inflection precludes an accurate test for any 
3.25 pm feature using our current method. The ice band optical depth plot (Fig. 5a) shows 
these crystalline ice features clearly. The derived optical depth is t(3.1) = 2.4 ± 0.2, higher 
than the previous estimate cited above, presumably because the earlier measurements were 
contaminated by light from W51/IRS2(W). The true optical depth toward the source might 
be even higher if scattered light or light from unresolved sources which does not undergo the 
same amount of ice extinction is present in the aperture. This would not be unreasonable 
since the source is quite distant (d=7500 pc, Goldader & Wynn- Williams 1994) and even a 
small effective aperture will subtend an extended area at the source. 



3.2. RNO 91 



RNO 91, in the L43 dark cloud, is a young stellar object which illuminates a reflection 
nebula. It has a 3 pm H 2 ice absorption band previously estimated to have an optical 
depth r — 1.29 from a 2.8-3.8 pm spectrum (Weintraub et al. 1994). The K(2.2 pm) band 
flux is extended in two lobes (Weintraub et al. 1994). All of our spectra were taken at the 
peak 3.5 pm position, presumably the exciting star. The K band spectrum (Fig. 4b) shows 
hot CO gas overtone absorption longward of 2.3 pm. 

To estimate the ice band depth, a 3rd order polynomial (in the log) was drawn between 
2.1-2.3 and 3.8-3.9 /im, avoiding the CO absorption (Fig. 4b). The derived ice band optical 
depth profile is shown in Fig. 5b. The H 2 ice toward this source is mostly amorphous. 
The derived optical depth is t(3.1) = 2.1 ± 0.2, higher than the estimate of Weintraub et 
al. (1994). Our optical depth is probably a better estimate as the 2 /im spectral continuum 
was included. 

This is the first published 2 /jm spectrum of this object to our knowledge. Because 
the goal was simply to define the continuum, the integration time was short and the 
signal-to-noise ratio low. However CO gas overtone absorptions can clearly be seen. 
Many young stellar objects show strong CO overtone absorption, comparable to late-type, 
main-sequence giants, though they have low total luminosities (Casali & Eiroa 1996; Greene 
& Lada 1996). The combined equivalent width of the CO 2-0 and 4-2 bands at 2.29 and 
2.35 /im in RNO 91, defined as in Greene & Lada (1996), is « 15 A. Two other photospheric 
features, the atomic Na and Ca lines at 2.21 and 2.26 /im, are weak, with equivalent widths 
^ 2 A. With these line strengths, RNO 91 falls between main sequence dwarf and giant 
photospheres, as do several young stellar objects in the p Oph cloud (Greene & Lada 



- 8- 



1996). As the bolometric luminosity of RNO 91 is low (« 3.7 L , Chen et al. 1995), the 
observations are consistent with photospheric absorption by a young star characterized 
by low surface gravity. Higher resolution spectra of this object should be able to confirm 
whether the CO absorption bands are in fact photospheric or due to a circumstellar disk. 



4. Determination of Spectral Features 



The method of defining the local continuum in order to extract the 3.25 and 3.47 
//m features is discussed in Sellgren, Smith, & Brooke (1994), Sellgren et al. (1995), and 
Brooke, Sellgren, & Smith (1996). To summarize the standard case, the local continuum 
was defined to be represented by points in the ranges 3.13-3.17 /mi, 3.33-3.37 /tm, and 
points longward of 3.61 fim. The continuua were estimated using low order polynomial fits 
to the logarithms of the flux densities of the continuum points, with the constraint that the 
best fit fall above those points held to contain features and look "reasonable" to the eye. In 
most cases, the best fit was a single 3rd order polynomial in the log or the joining of two 
3rd order polynomials. In cases where there were clearly no strong broad 3.47 or solid 
CH3OH features, the continuum was allowed to extend down to 3.58 to improve the 
baseline at the longer wavelengths. The adopted continuua are shown in Figs. 1-4. 

An exception to this standard case was made when the 3.53 jum absorption feature of 
solid CH3OH was clearly present: W33A, GL 2136, and MonR2/IRS2. The C-H stretch 
bands of solid CH3OH absorb strongly between approximately 3.33 /im and 3.65 /iin 
(Hudgins et al. 1993). For the sources with solid CH3OH, only points near 3.33 /im were 
used to define the continuum in that region. At the long wavelength end, our spectra lacked 
sufficient wavelength coverage. We used data from previously published spectra to help 
define the continuum, after first scaling the earlier spectra to the appropriate flux level. 
The earlier spectra were from Allamandola et al. (1992) for W33A, Schutte et al. (1996a) 
for GL 2136, and Smith, Sellgren & Tokunaga (1989) for MonR2/IRS2. The data taken 
from these spectra are indicated in Figs, la, lb, and 3a. The points used to define the 
continuum were from 3.65 /im to 3.70 fim. 

The resulting optical depths for sources with no detected CH 3 OH are shown in Figs. 
6 and 7. The 3.47 /im feature and the 3.25 /im feature (where detected) were fitted with 
gaussians with the central wavelength, peak optical depth, and full-width at half-maximum 
(FWHM) as free parameters. Uncertainties in these parameters were estimated from the 
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standard deviations of the results for fits using other baselines that looked reasonable. 
These are larger than the formal statistical uncertainties associated with the point-to-point 
scatter. Tables 2 and 3 summarize the results. There is no particular significance attached 
to the use of gaussians; the purpose of the fitting is only to extract reasonable estimates 
of the parameters of symmetric features. For W33A, GL2136, and MonR2/IRS2, the fits 
in the 3.47 fim region consisted of the sum of a gaussian and a CH 3 OH ice template, as 
discussed below. 

Although a local continuum could be drawn in the 3.25 //m region in most cases by 
our technique, the cancellation of telluric features was often insufficient to put any useful 
constraints on a 3.25 /im feature. One exception was Elias 29. No feature was detected but 
a 3a limit of r(3.25) < 0.03 was estimated (Fig. 7c). 

It is important to note that the actual absorption profiles of the species responsible 
for the features may extend further up into the long wavelength wing. Since we fit only a 
local continuum, our technique is sensitive only to the excess absorption at 3.47 /im or 3.25 
/im below this continuum, and the derived optical depths may be lower limits to the true 
contributions of the absorbers. 



5. Notes on Some Individual Sources 



W33A- Our estimate of the optical depth of the 3.47 /im feature is 0.29 ±0.04, considerably 
higher than the value of 0.15 estimated by Allamandola et al. (1992). This is due to the 
fact that the points we held to be continuum at the long-wavelength end were at longer 
wavelengths to take into account the presence of solid CH 3 OH. We continue to consider 
this value to be a lower limit to the true line-of-sight value, as the shape of the 3.1 /im 
H 2 ice band suggests that light suffering less ice extinction, possibly scattered light near 
the source or unresolved sources, is filling in the band to some extent (Brooke, Sellgren, & 
Smith 1996). The 3.25 /xm region was too faint to fit due to the deep H 2 ice band. 

GL 2136- The H 2 ice band optical depth toward this source has been recently measured 
to be r(3.1) = 3.2 ± 0.2 (Kastner & Weintraub 1996). The H 2 ice is partially crystalline 
(Fig. lb). The crystalline ice inflection at 3.20 /tm precludes an accurate test for the 3.25 
/im feature with our current technique. 

SVS20 SVS20 is a double source in the Serpens cloud, with a mostly N-S separation of 
1.6" (Eiroa et al. 1987). Our observations are of SVS20(S), which is 3 times brighter at 3.6 
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yum (Huard, Weintraub & Kastner 1997). Both objects exhibit H 2 ice absorption bands 
(Eiroa & Leinert 1987; Huard, Weintraub & Kastner 1997). 

RCrA/IRSl-Graham (1998) estimated a 3.47 fiva. feature optical depth of m 0.06. Our 
value is 0.041 ± 0.005, not remarkably different. 

RCrA/IRS2-Graham (1998) also estimated a 3.47 jj,m feature optical depth of ~ 0.06 for 
this source (TS 13.1). Our value is 0.041 ± 0.011, not remarkably different. 

GL 2591 - Although this source is relatively bright, our spectrum suffers from a poor 
cancellation of telluric features. No attempt was made to fit a continuum to the 3.25 fj,m 
region. 



6. The 3.47 fim Feature 



Using all of the available data, the best estimates of the central wavelength and width 
of this feature are A = 3.469 ± 0.002 fim (2883 cm" 1 ) and FWHM=0.105 ± 0.004 fim (87 
cm -1 ). Table 4 summarizes all of the detections of the feature to date. Figures 8a and 8b 
demonstrate that the 3.47 fim absorber is better correlated with the 3.1 fim H 2 ice band 
depth, than with the 9.7 fim silicate absorption band depth. References for the H 2 ice and 
silicate optical depths are given in Table 4. This is the same comparison done by Brooke, 
Sellgren, & Smith (1996), but the number of data points is now roughly doubled. The best 
weighted linear fit, exclusive of the upper limit, is 

r(3.47) = (0.033 ± 0.002) r(3.1) - (0.004 ± 0.004) (1) 

which passes near the origin. These coefficients are consistent with the values derived by 
Brooke, Sellgren, & Smith from the more limited data. The linear correlation coefficient 
is r = 0.88. The weighted fit is especially sensitive to the optical depth of GL961E (open 
circle at r(3.1) = 2.46 in Fig. 8b), as its error as estimated by Brooke, Sellgren, & Smith 
is significantly lower than the points around it. In the event that there may be some as 
yet undetermined systematic effects which affect the error determinations, we also give the 
unweighted fit: 

r(3.47) = 0.042 r(3.1) -0.010 (2) 

In any case, the abundance of the 3.47 fim absorber is closely related to that of H 2 ice 
over a wide range in ice extinction. 
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If in fact the 3.47 fim feature is due to the stretching vibration of sp 3 C-H bonds, 
the linear relation of Fig. 8b corresponds to roughly one C-H bond for every two H 2 
molecules, and a possible explanation for the correlation is that both H 2 ice and the C-H 
bonds form in step on grains by hydrogen-addition reactions (Brooke, Sellgren, & Smith 
1996). Detections of longer-wavelength C-H vibrational modes are needed to secure the 
identification; part of the broad 6.8 jum absorption feature may be due to as yet unidentified 
C-H deformation modes (Schutte et al. 1996b). 



7. The 3.25 Feature 



The best estimates of the central wavelength and width of this feature are 
A = 3.250 ± 0.004 /mi (3077 cm" 1 ) and FWHM=0.069 ± 0.006 pm (74 cm" 1 ). Figures 
9a and 9b show the optical depth of the 3.25 /j,m feature compared to silicates and H 2 
ice. With so few points, it is impossible to draw a firm conclusion, but the data to date, 
including the upper limit, suggest a better correlation with the silicates. (The best linear 
fit with the silicates has a nominal linear correlation coefficient r = 0.90, but with only 5 
points, there is a ~ 5 % probability that the data are completely uncorrelated.) 

If additional data show a correlation of the 3.25 /im feature with the silicate feature, 
this would provide evidence to support the Sellgren et al. (1995) identification of the 
feature as the C-H stretch band of aromatic hydrocarbons at low temperature. Aromatic 
hydrocarbons are thought to be ubiquitous in the interstellar medium and should correlate 
better with the refractory dust in molecular clouds than with volatile ices. 

The C-H stretch absorption band of aromatic hydrocarbons should be accompanied 
by absorption counterparts of the well-known aromatic emission features near 6.2, 7.7, 
8.6, and 11.3 [im. A possible absorption feature at 6.25 jum toward the young stellar 
object NGC7538/IRS9 has been suggested to be the C-C stretching mode of aromatic 
hydrocarbons (Schutte et al. 1996b). However, the identification and optical depth are 
still uncertain due to the presence of other absorbers. There are as yet no published 
high-sensitivity spectra of this object in the 3.25 ^m region. 

A test of the possible correlation of the 3.25 fxm feature with the refractory dust would 
be the clear detection of the feature along a line of sight with heavy dust extinction but 
little H 2 ice, as can be found in the diffuse interstellar medium (ISM). Such a detection 
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has not yet been made. Previous searches for the aromatic C-H stretch absorption in 
the diffuse ISM have been made near 3.29 /zm, the wavelength of the aromatic emission 
feature. Pendleton et al. (1994) discussed a possible absorption feature at 3.28 /im, perhaps 
due to aromatic hydrocarbons, toward three Galactic Center sources whose extinction is 
believed to be dominated by dust in the diffuse ISM. But the definition of the continuua 
was difficult, so they placed upper limits of r(3.28) <; 0.02 on this feature. Schutte et al. 
(1998) also placed an upper limit of r(3.3) ^ 0.02 toward one Wolf-Rayet star in their 
sample of diffuse ISM sources. Sellgren et al. (1995) proposed that the 3.25 /im feature 
observed in absorption in molecular clouds may be due to the same aromatic hydrocarbons 
responsible for the 3.29 /im emission feature, but shifted to a shorter wavelength due to a 
difference in temperature between cold absorbers and hot emitters. High-sensitivity spectra 
of diffuse ISM sources in the 3.25 /im region are needed to test this proposal. 

However an absorption feature at 6.2 /im identified as the aromatic C-C stretch has 
recently been identified toward 5 Wolf-Rayet stars and 2 Galactic Center sources in Infrared 
Space Obsevatory (ISO) spectra (Schutte et al. 1998). The extinction to all of these objects 
is believed to come primarily from the diffuse ISM. The optical depth of the 6.2 /im feature 
is well-correlated with the 9.7 /im silicate absorption optical depth toward these sources. 
In Figure 10, we compare the relationships of the possible aromatic absorptions at 6.2 and 
3.25 /im and silicate absorption. For this comparison, we multiplied the r(6.2) values of 
Schutte et al. by 0.25, the factor needed to bring the best-fit linear slopes into agreement. 
Figure 10 suggests that the apparent relationship between r(3.25) and r(9.7) is consistent 
with the linear correlation of r(6.2) vs. r(9.7), after correcting by a scale factor. This 
agreement is what one would expect if the abundance of aromatic hydrocarbons relative to 
silicates, and the relative strengths of C-H stretch and C-C stretch modes, were similar in 
both molecular clouds and the diffuse ISM. 

Assuming that the 3.25 /im feature is due to aromatic hydrocarbons, abundances 
can be estimated following Sellgren et al. (1995). The abundances of aromatic C-H 

bonds relative to hydrogen gas lie in the range N(C-H)/N(H)« 1 3 x 10~ 5 (see Table 

5). Adopting a value for the ratio of aromatic carbon atoms in hydrogen bonds to total 
aromatic carbon atoms of / = 0.4 (Sellgren et al. 1995), roughly 8-23 % of cosmic carbon 
would be in aromatic hydrocarbons. 



8. Solid CH 3 OH 
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The 3.53 /tin feature due to solid CH 3 OH was detected toward three objects here: 
W33A, GL2136, and MonR2/IRS2. The detection toward MonR2/IRS2 is new. A previous 
spectrum of NGC7538/IRS9 showing solid CH 3 OH (Brooke, Sellgren, & Smith 1996) is also 
analyzed in detail here. For consistency with the other sources, the optical depth profile 
of NGC7538/IRS9 in the 3.47 /zm region was re-derived using the same definition of the 
local continuum. This required some additional data at the long-wavelength end; these 
points were obtained from Allamandola et al. (1992) and were corrected for a slight slope 
difference. The derived gaussian parameters are given in Table 2. They are not significantly 
different from those derived earlier. We thus present analyses of four sources in which the 
3.53 /iin feature of solid CH 3 OH has been firmly identified to date. The feature has also 
recently been detected toward a fifth source, RAFGL7009S (Dartois et al. 1998). 

The peak wavelength of the CH 3 OH z/ 3 band at 3.53 /im can be used to constrain the 
composition of the methanol-containing ices (Schutte et al. 1996a). In mixtures with H 2 
ice, the feature shifts to shorter wavelengths (higher frequency) with increasing H2O fraction. 
We used laboratory spectra of pure CH 3 OH and a 100:10:1:1 H 2 0:CH 3 OH:CO:NH 3 mix 
("Weak Interstellar Mix"), both at 10K, obtained in digital form from the AAS CD-ROM 
version of Hudgins et al. (1993). The peak in pure CH 3 OH occurs at 3.535 /im (2829 cm -1 ). 
In the weak ISM mix, it occurs at 3.530 /im (2833 cm -1 ). 

In order to fit the spectra, local continuua similar to those used for the observations 
were applied to the CH 3 OH laboratory data. Each of the spectra in Fig. 11 was 
simultaneously fitted with a gaussian for the 3.47 /zm feature and a CH 3 OH ice template. 
Optical depths are quoted for the 3.53 /im feature alone. An additional narrow absorption 
feature near 3.47 /zm in W33A discussed below was deleted from the fit for this object. 

The features in GL2136 and NGC7538/IRS9 occur at 3.536 /im (2828 cm" 1 ), in good 
agreement with pure CH 3 OH ice. In W33A, the feature occurs at 3.534 /xm (2830 cm -1 ). 
We have fit this feature with a mix of the two samples described above, as indicated in 
Table 6. 

The feature in MonR2/IRS2 occurs at 3.528 /im (2834 cm" 1 ). This is too short a 
wavelength to be consistent with pure CH 3 OH. The weak ISM mix provides a better, though 
still not exact, fit with r(3.53) = 0.033 ± 0.010. Presumably an even more water-rich mix 
would improve the fit. It would be desirable to have laboratory spectra of such mixtures. 

Previous studies indicated that the solid methanol toward W33A, GL 2136, and 
NGC7538/IRS9 is in a CH 3 OH-rich phase with CH 3 OH/H 2 £ 0.5 (Skinner et al. 1992; 
Allamandola et al. 1992; Schutte et al. 1996a; Dartois et al. 1998) and our fits confirm 
this. The 3.53 /im feature in MonR2/IRS2 suggests that solid methanol can also be found 
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within H 2 0-rich ice having CH 3 OH/H 2 0^ 0.1. 

An additional indicator of CH 3 OH in H 2 0-rich environments would be the presence of 
two weak features which appear at 3.38 yum and 3.41 fxm (Schutte et al. 1996a). These can 
be seen in the laboratory data in Fig. 11. The present signal-to-noise is not high enough to 
constrain the presence of these features in any of the sources. 

Column densities and CH 3 OH abundances are listed in Table 6. For W33A, the column 
density of H 2 ice is uncertain. The range given results from converting the observed 
optical depths of the 3.1 and 6.0 yum features to column densities. But, as discussed above, 
the 3.1 fxm band is most likely filled in by light which undergoes less ice extinction. On 
the other hand, the 6.0 yum band may contain other absorbers (Schutte et al. 1996b). The 
range given should bracket the true line-of-sight column density. It should be noted that 
the CH 3 OH column density may be an underestimate if the filling-in noted above is also 
important at 3.53 yum. 

In the other three sources, the CH 3 OH/H 2 ice abundance derived from the 3.53 fj,m 
band ranges from 3 to 5 %. The derived value of the abundance for GL 2136 is in agreement 
with that derived by Schutte et al. (1996a) from the same band. The derived values for 
W33A and NGC7538/IRS9 are roughly a factor 3 lower than those of Allamandola et al. 
(1992), due to the use of different band strengths, different adopted continuua, and different 
assumed contributions of the 3.47 /im. feature. Our abundance range for W33A, 3-16 % is 
in reasonable agreement with that of Dartois et al. (1998), 5-22 %. 

Upper limits to the CH 3 OH/H 2 ice ratio toward other sources have been estimated 
(Chiar, Adamson, & Whittet 1996; Dartois et al. 1998). The most stringent is that for 
W3/IRS5 for which the upper limit is of order 1 %. In contrast, Dartois et al. estimate 
a very high abundance of 30 % for RAFGL7009S. Thus solid CH 3 OH is common along 
molecular cloud lines-of-sight with H 2 ice, with a typical abundance of a few percent, but 
its abundance definitely varies. 

There are two notable aspects about the methanol abundances. First, the typical 
values (1-5 % relative to H 2 0) of the CH 3 OH abundance are very similar to the values 
in comets, which are derived from the same band in gas-phase fluorescence (Hoban et al. 
1991; Bockelee-Morvan, Brooke, & Crovisier 1995; DiSanti et al. 1995). This provides some 
supporting evidence for the direct incorporation of interstellar ice grains into comets. 

Second, the CH 3 OH/H 2 ice ratio toward NGC7538/IRS9 is similar to the ratios 
toward GL 2136 and MonR2/IRS2, but the relative abundance of solid CO trapped in 
non-polar ices toward this source is ;> 5 times higher (Tielens et al. 1991; Chiar et al. 
1998). Thus the abundance of solid CH 3 OH is not tightly coupled to this CO component, 
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which is believed to form in dense regions of molecular clouds where most of the hydrogen is 
molecular (Tielens et al. 1991). This favors theories which postulate the formation of solid 
CH3OH by hydrogen-addition reactions in regions where atomic hydrogen predominates, as 
has also been proposed for H 2 ice (e.g. d'Hendecourt, Allamandola, & Greenberg 1985). 
The abundance of solid CO trapped in polar ices (presumably H2O) is also similar in all 
three sources (2-3 %, Tielens et al. 1991) and this CO component may form under similar 
conditions. 



9. Additional Narrow Absorption Near 3.47 /im 



There is an additional narrow absorption feature near 3.47 /im (2882 cm -1 ) in the 
spectrum of W33A (Fig. 11). This region was excluded from the fit with the 3.47 /im 
gaussian plus solid methanol for W33A. It is unlikely to be simply structure within the 
broad 3.47 //m feature as it does not appear at the same level in well-measured 3.47 //m 
features like those toward MonR2/IRS3 (Sellgren et al. 1995) and W3/IRS5 (Brooke, 
Sellgren, & Smith 1996). The feature is detected at the 3a level. A similar absorption 
feature was noted by Dartois et al. (1998). 

Before examining the feature toward W33A, we discuss the spectrum of GL 2136 in 
this region. Schutte et al. (1996a) proposed that a feature at 3.473 /im (2879 cm -1 ) with 
FWHMsi 25 cm -1 was present in GL 2136 with an optical depth r « 0.01. They identified 
the feature as the z/ 5 band of solid H 2 CO in a mixture with H 2 and CH 3 OH. Our spectrum 
of GL 2136 does not show the same structure which led Schutte et al. to claim the presence 
of a feature, and the current spectrum has higher signal-to-noise. In our spectrum, this 
region resembles the peak of a broad 3.47 //m absorption feature like that seen in other 
sources. However due to uncertainties in baselines we cannot rule out some contribution 
from solid H2CO at the bottom of the band. Our 3cr upper limit is r <^ 0.01, i.e. comparable 
to the claimed detection, which corresponds to a solid H 2 CO column density of ^ 3 x 10 17 
cm -2 and an abundance of ^ 6 % relative to H 2 using the same band width and strength 
adopted by Schutte et al. If this much H 2 CO were present, part of the 3.53 /im band would 
be due to the v\ H 2 CO band and the solid CH 3 OH abundance would be roughly 9 % lower. 
But as we see no evidence for solid H 2 CO, we prefer attributing all of the 3.53 /im band to 
solid CH3OH. 

The narrow 3.47 /im feature in W33A falls at a wavelength consistent with solid H 2 CO 
in various mixtures with CH 3 OH and H 2 (Schutte et al. 1996a), but the feature appears 
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asymmetric, unlike H 2 CO. If the excess absorption is interpreted as entirely due to the 
u 5 band of H 2 CO, the peak optical depth of r ~ 0.06 would imply a column density of 
1.8 x 10 18 cm -2 . As discussed above, the H 2 ice column density to W33A is uncertain, 
but this upper limit to the optical depth corresponds to a solid H 2 CO abundance of ^4-20 
% relative to H 2 0. If this feature were due to H 2 CO, the v\ H 2 CO band would contribute 
to the 3.53 {xm band and methanol abundances would be ~ 12 % lower. 

Another possible contributor to the narrow 3.47 /im feature is solid ethane, C 2 H6, 
which has a narrow (FWHM^ 6 cm -1 ) band at 3.472 /im (2880 cm -1 ) (Boudin, Schutte, & 
Greenberg 1998). A stronger C 2 H 6 band occurs at 3.365 /xm (2972 cm -1 ) and the optical 
depth of this feature should be 2-3 times higher. The present data for W33A would be 
consistent with a 3a optical depth upper limit of approximately 0.06 for this feature. Thus 
up to about half of the apparent narrow 3.47 ^m feature could be due to solid C 2 H 6 . This 
would correspond to a solid C 2 H 6 column density of 3.0 x 10 16 cm -2 or 0.1-0.3 % relative to 
n 2 ice. This is near the abundance of C 2 H 6 in comets C/1996 B2 Hyakutake and C/1995 
01 Hale-Bopp of ~0.4 % relative to H 2 (Mumma et al. 1996; Weaver et al. 1998). Future, 
higher signal-to-noise spectra of W33A in this region should be able to distinguish possible 
contributions of solid H 2 CO and C 2 H 6 . 



10. Conclusions 

1) The optical depth of the broad absorption feature at 3.47 fim seen on the long-wavelength 
wing of molecular cloud H 2 ice bands is correlated with the optical depth of H 2 ice, in 
agreement with an earlier conclusion based on a smaller data set. If the feature is due to 
C-H bonds, the abundance of the C-H bonds must be closely related to that of H 2 ice. 

2) The broad absorption feature at 3.25 ^m was detected toward two more young stellar 
objects. There are insufficient data to draw any firm conclusions, but the optical depths 
measured to date correlate better with the refractory silicates than with H 2 ice, providing 
evidence that supports an identification with aromatic hydrocarbons at low temperature. 

3) The 3.53 /mi solid CH 3 OH features toward W33A, GL2136, and NGC7538/IRS9 can be 
fit by CH 3 OH-rich ices, while the wavelength of the feature toward MonR2/IRS2 suggests 
an H 2 0-rich ice environment. Solid methanol abundances toward GL 2136, NGC7538/IRS9, 
and MonR2/IRS2 are 3-5 % relative to H 2 ice, similar to cometary abundances. 

4) There is an additional narrow absorption feature near 3.47 fim toward W33A. Possible 
contributors are solid H 2 CO and solid C 2 H 6 . 
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5) The H 2 ice toward W51/IRS2 is located primarily in front of the eastern component, is 
partially crystalline, and has an optical depth r(3.1) = 2.4 ± 0.2. Spectra of RNO 91 from 
2 to 4 /im reveal stellar (or circumstellar) CO gas absorption and an H 2 ice band optical 
depth r(3.1) = 2.1 ±0.2. 
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Notes to Table 1: 
a Integration time 

b Assumed blackbody temperature of standard star. 
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Table 1: Log of Observations 



Object 


UT 


Tel 


Range 


t a 


Standard 


Spectral 


u 

rri O 

Jbb 










(H 


(sec) 


Star 


Type 


(K) 


MonR2/IRS2 


12 Jan 


1996 


2.1m 


3.3-3.6 


180 


e Ori 


B0I 


26000 


p Oph/Elias 29 


12 Jun 


1996 


UKIRT 


3.1-3.4 


288 


49 Lib 


F8V 


6200 




13 Jun 


1996 




3.3-3.6 


288 


BS 5923 


F6V 


6400 




16 Jul 


1996 




2.9-3.2 


288 


5 Sco 


B0IV 


27000 


p Oph/WL 6 


13 Jun 


1996 


UKIRT 


3.1-3.4 


432 


BS 6496 


F7V 


6300 




13 Jun 


1996 




3.3-3.6 


432 


BS 6310 


F3V 


6700 




16 Jul 


1996 




2.9-3.2 


432 


5 Sco 


B0IV 


27000 


p Oph/Elias 21 


12 Jun 


1996 


UKIRT 


3.1-3.4 


288 


BS 6310 


F3V 


6700 




12 Jun 


1996 




3.3-3.6 


288 


BS 6310 


F3V 


6700 




16 Jul 


1996 




2.9-3.2 


432 


5 Sco 


B0IV 


27000 


p Oph/Elias 33 


12 Jun 


1996 


UKIRT 


3.1-3.4 


288 


BS 6310 


F3V 


6700 




12 Jun 


1996 




3.3-3.6 


288 


BS 6310 


F3V 


6700 




16 Jul 


1996 




2.9-3.2 


432 


BS 7152 


F0V 


7200 


RNO 91 


06 Jun 


1996 


2.1m 


2.0-2.4 


480 


v Ser 


A2V 


9000 




06 Jun 


1996 




3.0-3.9 


240 


v Ser 


A2V 


9000 




13 Jun 


1996 


UKIRT 


3.1-3.4 


432 


BS 6496 


F7V 


6300 




13 Jun 


1996 




3.3-3.6 


432 


BS 6012 


F3V 


6700 




16 Jul 


1996 




2.9-3.2 


432 


5 Sco 


B0IV 


27000 


W33A 


06 Sep 


1995 


UKIRT 


2.9-3.6 


1440 


BS 6496 


F7V 


6300 


GL 2136 


06 Sep 


1995 


UKIRT 


2.9-3.6 


576 


BS 7126 


F5V 


6440 


Ser/SVS20 


04 Jun 


1996 


2.1m 


3.3-3.6 


180 


S Aql 


F2IV 


6890 




16 Jul 


1996 


UKIRT 


2.9-3.2 


240 


BS 6797 


F5V 


6440 




16 Jul 


1996 




3.1-3.4 


432 


BS 6797 


F5V 


6440 


RCrA/IRSl 


12 Jun 


1996 


UKIRT 


3.3-3.6 


288 


BS 7152 


F0V 


7200 




13 Jun 


1996 




3.1-3.4 


432 


BS 7152 


F0V 


7200 




16 Jul 


1996 




2.9-3.2 


288 


BS 7152 


F0V 


7200 


RCrA/IRS2 


13 Jun 


1996 


UKIRT 


3.1-3.4 


432 


BS 7152 


F0V 


7200 




13 Jun 


1996 




3.3-3.6 


432 


BS 7152 


F0V 


7200 




16 Jul 


1996 




2.9-3.2 


240 


5 Sco 


B0IV 


27000 


RCrA/IRS5 


13 Jun 


1996 


UKIRT 


3.3-3.6 


432 


BS 7152 


F0V 


7200 


W51/IRS2 


06 Jun 


1996 


2.1m 


2.0-2.4 


120 


15 Vul 


A7m 


7850 




06 Jun 


1996 




3.0-3.9 


320 


111 Her 


A3V 


8720 




12 Jun 


1996 


UKIRT 


3.1-3.4 


288 


BS 7172 


F8V 


6200 




13 Jun 


1996 




2.9-3.2 


432 


BS 7550 


F5V 


6440 


GL 2591 


06 Sep 


1995 


UKIRT 


2.9-3.6 


576 


r Cyg 


F3IV 


6700 
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Notes to Table 2: 

a Central wavelengths, full widths at half maximum, and peak optical depths of 
absorption features from gaussian fits. 

b Entries in parentheses are la uncertainties obtained from standard deviations of 
results using several different baselines (see text). 

c Interpreted as a lower limit (see text). 



- 21 - 



Table 2: 3.47 pm Absorption Feature Parameters 



Object 


Ao 


AA 


r 




(H 


(pm) 




MonKz/lKoz 


3.460 


n i no 

0.103 


0.083 




(0.006) 


(0.007) 


(0.013) 


p Oph/Ehas 29 


o a aa 

3.490 


0.113 


A A OA 

0.089 




(0.006) 


(0.009) 


(0.007) 


p Oph/WL o 


O A TO 

3.472 


0.111 


A A/^O 

0.062 




(0.003) 


(0.009) 


(0.012) 


p Oph/Ehas 21 


O AHA 

3.464 


a i a a 

0.104 


A A A O 

0.042 




(0.003) 


(0.012) 


(0.012) 


p Oph/Ehas 33 


3.457 


r\ i an 

0.120 


A A A /I 

0.024 




(0.009) 


(0.007) 


(0.006) 


T") \T/ \ f\-t 

KINO 91 


O A PC\ 

3.469 


0.118 


A A o r 

0.085 




(0.003) 


/A A "1 o\ 

(0.012) 


(0.015) 


TT TO O A 

W33A 


3.477 


A AAA 

0.099 


A OA AT 

0.290 




/a a i a \ 

(0.010) 


/A A "1 o\ 

(0.012) 


/a a a a \ 

(0.040) 


OL 2136 


3.471 


0.118 


A 1 QT 

0.137 




(0.010) 


/ A A A r - \ 

(0.005) 


/ A AO A \ 

(0.020) 


ber/b Vb20 


3.463 


0.111 


A AO A 

0.030 




(0.011) 


(0.009) 


(0.009) 


RCrA/lRbl 


3.469 


A 1 OO 

0.122 


A A A 1 

0.041 




(0.012) 


( r\ ai r*\ 

(0.016) 


/a a a r \ 

(0.005) 




^ 479 


D DQ9 


U.U11 




(0.005) 


(0.005) 


(0.011) 


RCrA/IRS5 


3.470 


0.111 


0.090 




(0.006) 


(0.005) 


(0.019) 


W51/IRS2(E) 


3.471 


0.094 


0.143 




(0.003) 


(0.005) 


(0.012) 


GL 2591 


3.464 


0.131 


0.045 




(0.005) 


(0.007) 


(0.005) 


NGC7538/IRS9 


3.485 


0.118 


0.130 




(0.003) 


(0.024) 


(0.017) 
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Notes to Table 3: 

a Central wavelengths, full widths at half maximum, and peak optical depths of 
absorption features from gaussian fits. 

b Entries in parentheses are la uncertainties obtained from standard deviations of 
results using several different baselines (see text). 

c 3a upper limit. 



Table 3: 3.25 /im Absorption Feature Parameters 



Object 


A 


AA 


r 




(/im) 


(/im) 




Ser/SVS20 


3.245 


0.061 


0.021 




(0.003) 


(0.007) 


(0.003) 


RCrA/IRSl 


3.239 


0.059 


0.032 




(0.003) 


(0.012) 


(0.010) 


p Oph/Elias 29 






<0.03 c 
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Notes to Table 4: The optical depths given in the table are for silicates at 9.7 pm; 
H 2 ice at 3.1 /im; and broad features at 3.47 and 3.25 yum. The silicate optical depth for 
W51/IRS2(E) was corrected for intrinsic silicate emission following Willner et al. 1982. 

a Willner et al. 1982. 

b Smith, Sellgren, & Tokunaga 1989. 

c Brooke, Sellgren & Smith 1996. 

d Whittet et al. 1988. 

e Hanner, Brooke, & Tokunaga 1998. 

f Sato et al. 1990. 

9 Chiar, Adamson, & Whittet 1996. 

h Smith, Sellgren, & Brooke 1993. 

* This work. 

j Sellgren, Smith & Brooke 1994. 

k Hanner, Brooke, & Tokunaga 1995. 

' Tanaka et al. 1990. 

m Kastner & Weintraub 1996. 

n J. Guertler, priv. comm. 

Eiroa & Hodapp 1989. 

p Whittet et al. 1996. 

q Tanaka et al. 1994. 

r Genzel et al. 1982. 
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Table 4: Summary of Absorption Feature Optical Depths 



Object 


r(9.7) 


r(3.1) 


r(3.47) 


r(3.25) 


W3/IRS5 


7.64 a 


3.48 6 


0.13 c 


— 


Tau/Elias 18 


0.43 d 


0.80 d 


0.022 c 


— 


HL Tau 


0.56 e 


0.77 f 


0.035 9 


- 


Tau/Elias 16 


0.66 d 


1.26 h 


0.032^ 


- 


BN 


3.28 a 


1.78 b 


0.034 c 


- 


Mon R2/IRS2 


— 


2.54 6 


0.083* 


— 


Mon R2/IRS3 


4.30 a 


1.14 6 


0.036 J ' 


0.049 J ' 


S255/IRS1 


5.11 a 


1.48 6 


0.024 c 


- 


GL961E 


2.11 a 


2.46 6 


0.068 c 


— 


p Oph/Elias 29 


1.51 fe 


1.80' 


0.089*' 


<0.03* 


p Oph/WL 6 


1.22 fe 


2.1' 


0.062* 


- 


p Oph/Elias 21 


0.61 fe 


0.77' 


0.042* 


— 


p Oph/Elias 33 


0.36 fe 


1.00' 


0.024*' 


— 


RNO 91 




2.F 


0.085* 




W33A 


7.84 a 


>5.4 a 


>0.29* 




GL 2136 


5.07 a 


3.2 m 


0.14* 




Ser/SVS20 


2.0 n 


1.0° 


0.030* 


0.021* 


RCr A/IRS 1 


1.2F 


1.44" 


0.041* 


0.032* 


RCrA/IRS2 




1.1« 


0.041* 




RCrA/IRS5 




2.3 q 


0.090* 




W51/IRS2(E) 


3.7 r 


2.4* 


0.14*' 




GL 2591 


4.14 a 


0.92 6 


0.045* 




S140/IRS1 


3.97 a 


1.28 6 


0.027 c 


0.036 c 


NGC 7538/IRS1 


6.38 a 


1.29 6 


0.052 c 


0.078 c 


NGC 7538/IRS 9 


4.46 a 


3.28 b 


0.13*' 
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Notes to Table 5: 

a Optical depth of 3.25 fim absorption feature, la error in parentheses. References for 
r(3.25) are in Table 4. 

6 Column density of aromatic C-H bonds calculated from N(C-H)=r(3.25)Az//A with 
A— 1.7 x 10~ 18 cm/C-H bond (Sellgren et al. 1995) and Au = 74 cm -1 , la error in 
parentheses. 

c Column density of hydrogen gas estimated from N(H)=3.5 x 10 22 x r(9.7) cm -2 
(Tielens et al. 1991). References for r(9.7) are in Table 4. 

d Abundance of aromatic C-H bonds relative to hydrogen gas. la error in parentheses. 
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Table 5: Summary of 3.25 yum Features 



r(3.25) a N(C-H) 6 N(H) C [C-H/H] d 

(10 18 cm' 2 ) (10 23 cm- 2 ) 



Mon R2/IRS3 


0.049 


2.1 


1.5 


1.4 




(0.007) 


(0.3) 




(0.2) 


Elias 29 


< 0.03 


< 1.3 


0.53 


< 2.4 


SVS 20 


0.021 


0.91 


0.70 


1.3 




(0.003) 


(0.13) 




(0.2) 


RCrA/IRSl 


0.032 


1.4 


0.42 


3.3 




(0.010) 


(0.4) 




(0.9) 


S140/IRS1 


0.036 


1.6 


1.4 


1.1 




(0.007) 


(0.3) 




(0.2) 


NGC7538/IRS1 


0.078 


3.4 


2.2 


1.5 




(0.013) 


(0.6) 




(0.3) 
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Notes to Table 6: 

a Optical depth of 3.53 /xm absorption feature, la error in parentheses. 

b For the Weak ISM Mix, a FWHM Au=29 cm' 1 and absorbance A=7.5 x KT 18 
cm/molecule (Hudgins et al. 1993) were assumed. For pure CH 3 OH, a FWHM Az^=29 
cm -1 and absorbance A=6.6 x 1CT 18 cm/molecule (Schutte et al. 1996a) were assumed. 

c Column density of solid CH 3 OH calculated from N(CH 3 OH)=r(3.53)Az//A 
(Allamandola et al. 1992). la error in parentheses. 

d Column density of H 2 ice. 

e Abundance of solid CH 3 OH relative to H 2 ice. la error in parentheses. 

f Estimated from the 3.1 /xm band measured by Smith, Sellgren, & Tokunaga (1989). 

9 Estimated from the 3.1 and 6.0 fxm bands by Allamandola et al. (1992). 

h Estimated from the 3.1 /im band by Schutte et al. (1996a). 

1 Estimated from the 3.1 /im band by Allamandola et al. (1992). 
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Table 6: Solid CH 3 OH Detections at 3.53 /im 





r(3.53) a 


Fit Type 6 


N(CH 3 OH) c 
(10 17 cm" 2 ) 


N(H 2 0) d 
(10 19 cm- 2 ) 


[CH 3 OH/H 2 0] e 


Mon R2/IRS2 


0.033 


Weak ISM Mix 


1.3 


0.42 / 


0.031 




(0.010) 




(0.4) 




(0.009) 


W33A 


0.27 


Pure CH 3 OH 


11.9 


0.9-4.2^ 


0.03-0.16 




(0.05) 


+ 


(2.2) 








0.07 


Weak ISM Mix 


2.7 








(0.02) 




(0.8) 
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Figure Captions 

Figure l(ab) — Spectra from 2.9 to 3.6 fim with resolution A/AA « 1200 (solid lines). 
Dashed lines are polynomial fits in local continuum regions (see text), (a) Filled circles are 
data from Allamandola et al. (1992) scaled to the same flux level, (b) Filled circles are data 
from Schutte et al. (1996a) scaled to the same flux level. Apparent hydrogen line emission 
features at 3.04 yum (HI Pfe) and 3.30 jj,m (HI PfS) in GL 2136 may be due all or in part to 
corresponding absorption features in the standard stars, which were not corrected for. 

Figure 2(a-h) — Spectra from 2.9 to 3.6 ^m with effective resolution A/AA ~ 1300 
(1200 for GL 2591) are shown as solid lines. Dashed lines are polynomial fits in local 
continuum regions (see text). Apparent hydrogen line emission features at 3.04 /xm (HI 
Pfe) and 3.30 /jm (HI Pf<5) in Elias 29, SVS20, RCrA/IRSl, and RCrA/IRS2 may be due 
all or in part to corresponding absorption features in the standard stars, which were not 
corrected for. 

Figure 3(ab) — Same as Figure 2 for objects observed only between 3.3 and 3.6 fxm. 
Filled circles for MonR2/IRS2 are data from Smith, Sellgren, & Tokunaga (1989) scaled to 
the same flux level. 

Figure 4(ab) — 2-4 /j,m spectra, (a) For W51/IRS2(E), spectral resolution was ~ 700 
in the 2.2 jj,m region, ps 1300 in the 2.9-3.4 /im region and « 530 in the 3.4-3.9 
region. Long dashed line is estimated continuum. A scaled blackbody (dashed-dotted line) 
indicates the approximate color temperature. Short dashed line is the local continuum 
for the 3.47 fim feature. Wavelengths of the indicated lines are (in microns): 2.059 (Hel 
2p 1 P°-2s 1 S), 2.113 (Hel 4s 3 S-3p 3 P°), 2.122 (H 2 v=l-0 S(l)), 2.166 (HI Br 7 ), 3.039 (HI Pfe), 
3.297 (HI Pf<5), 3.741(HI Pfy). After correction for hydrogen absorption in the standard 
stars, the line flux ratios (Wm -2 ) relative to Btj of the indicated lines are (in increasing 
wavelength order): 0.63, 0.04, 0.03, 1.0, 0.38, 0.91, 5.21. (6) For RNO 91 the effective 
spectral resolution was ~ 400 in the 2.2 jj,m region, 1300 in the 2.9-3.6 /im region, and 
200 in the 3.6-3.9 fim. region. Long dashed line is estimated continuum. Short dashed 
line is the local continuum for the 3.47 /im feature. The wavelengths of the CO band heads 
are (in microns): 2.293 (2-0), 2.323 (3-1), 2.352 (4-2), 2.383 (5-3). 

Figure 5(a6) — Derived H2O ice band profiles from the continuua indicated in Figure 
4. Wavelengths of the hydrogen emission lines are given in Fig. 4. 

Figure 6(a-h) — Derived 3.47 /im feature optical depth profiles from continuua 
indicated in Figs. 2-3. Dashed lines are gaussian fits. A representative ±1<t error bar is 
offset lower right. 
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Figure 7(o-c) — Derived 3.25 /jm and 3.47 jim feature optical depth profiles from 
continuua indicated in Fig. 2. Dashed lines are gaussian fits. Some ±1<t error bars are 
indicated. 

Figure 8(ab) — Optical depths for the 3.47 /im feature vs. (o) silicates and (b) H 2 
ice. Solid line is the best weighted linear fit. Dashed line is best unweighted linear fit. Open 
circles are from Sellgren, Smith, & Brooke (1994) and Brooke, Sellgren, & Smith (1996); 
X's from Chiar, Adamson, & Whittet (1996); filled circles from this work. 

Figure 9(ab) — Optical depths for the 3.25 /an feature vs. (a) silicates and (b) H 2 
ice. Open circles are from Sellgren, Smith, & Brooke (1994) and Brooke, Sellgren, & Smith 
(1996); filled circles from this work. 

Figure 10 — Optical depths of the 3.25 fim absorption feature measured toward 
molecular cloud sources, and scaled values for the 6.2 /xm absorption feature measured 
toward sources whose extinction is believed to come primarily from the diffuse ISM (Schutte 
et al. 1998), compared to silicate absorption optical depths. Open circles are from Sellgren, 
Smith, & Brooke (1994) and Brooke, Sellgren, & Smith (1996); filled circles from this work. 
The Schutte et al. sources are 5 Wolf-Rayet stars (WR118, WR112, WR104, WR98a, 
WR48a) and an upper limit for Cyg OB2 No. 12 (small filled squares); and 2 Galactic 
Center sources, GCS3 and GCS4 (small filled diamonds). The scaling factor of 0.25 brings 
the best-fit linear slopes into agreement. 

Figure 11 — Derived optical depths from continuum fits for sources with the 3.53 /jhi 
feature. Short dashed line is the best fit consisting of the sum of a 3.47 /jul feature gaussian 
(long dashes) and solid CH 3 OH (dotted-dashed) as described in text. Some representative 
±lcr error bars are indicated. 
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